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Abstract. Pure singlets are typically disfavored as dark matter candidates, since they gener- 
ically have a thermal relic abundance larger than the observed value. In this paper, we 
propose a new dark matter mechanism called "assimilation" , which takes advantage of the 
baryon asymmetry of the universe to generate the correct relic abundance of singlet dark 
matter. Through assimilation, dark matter itself is efficiently destroyed, but dark matter 
number is stored in new quasi-stable heavy states which carry the baryon asymmetry. The 
subsequent annihilation and late-time decay of these heavy states yields (symmetric) dark 
matter as well as (asymmetric) standard model baryons. We study in detail the case of pure 
bino dark matter by augmenting the minimal supersymmetric standard model with vector- 
like chiral multiplets. In the parameter range where this mechanism is effective, the LHC can 
discover long-lived charged particles which were responsible for assimilating dark matter. 
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1 Introduction 

The existence of dark matter (DM) is striking evidence for physics beyond the standard 
model (SM) [1-5]. Stable particles with Fermi-scale masses and electroweak annihilation 
cross sections, known as Weakly Interacting Massive Particles (WIMPs), are prime DM 
candidates since they have the right thermal relic abundance to account for observations 
[6-10]. Supersymmetry (SUSY) with conserved ii-parity naturally provides WIMP DM if 
the lightest super symmetric particle (LSP) is a neutralino [11-15]. This can occur in the 
Minimal Supersymmetric Standard Model (MSSM) as well as in various generalizations. 

Depending on the soft mass spectrum, it is often the case that the LSP is a nearly pure 
bino eigenstate. The bino is a gauge singlet, and its annihilation through p-wave-suppressed 
t-channel sfermion exchange is too small to yield the correct DM relic density. For this reason, 
a pure bino is typically excluded by overclosure, and mixed neutralino DM is often assumed 
in the MSSM. There are of course ways to accommodate bino DM. The bino might be almost 
degenerate in mass with the lightest sfermion such that thermal production can be dominated 
by co-annihilation [16-18]. Alternatively, one can choose a reheat temperature lower than the 
freezeout temperature to achieve the correct bino relic density [19]. Besides the MSSM bino, 
SM gauge singlets are well-motivated DM candidates, and the correct thermal abundance can 
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arise if the singlet interacts with portal fields that provide additional annihilation channels 
[20-24] , or if the singlet mixes with "higgsino-like" states [25-28] . 

In this paper, we propose an alternative mechanism to achieve the correct thermal abun- 
dance of singlet DM, by taking advantage of matter/antimatter asymmetries. We introduce 
a new DM process called "assimilation", which takes the schematic form 

^^cp. (1.1) 

Here, x is singlet DM, ^ and ^ are new quasi-stable heavy states which carry an asymmetry, 
and is a SM field. The tildes in x and ^ denote a multiplicatively conserved DM number, 
such as ii-parity in SUSY theories. In the presence of a asymmetry, the process in Eq. (1.1) 
is very effective at destroying singlet particles x- process, DM number (e.g. i?-parity) 

is conserved and stored in In this way, DM is "assimilated" into the heavier state ^. 

The reason an asymmetry is necessary is that ^ has to be sufficiently dense at tem- 
peratures of the order of the DM mass for Eq. (1.1) to be relevant. Without an asymme- 
try, ^ and ^ are quickly depleted by the crossed version of Eq. (1.1). Motivated by the 
baryon/antibaryon asymmetry, we will focus on the case where ^ and ^ are charged under 
[/(I) B and carry the full baryon asymmetry at early times. We will not give an explicit model 
of ^'-ogenesis in this paper, though certain standard baryogenesis mechanisms can be adapted 
to generate a ^ asymmetry. Note that unlike asymmetric DM [29-40], DM itself does not 
have an asymmetry. Unlike in co-annihilation [16], ^ and ^ are sufficiently long-lived that 
they can be considered stable during DM freezeout.^ 

Since ^ are not gauge singlets, they typically have large annihilation cross sections, 
especially in SUSY theories. Assuming that ^ is heavier than ^, the process 

^ ^r^' (1.2) 

efficiently depletes DM number until ^ freezes out. Analogous to the super- WIMP scenario 
[42, 43], at late times, the decay of ^' produces singlet particles x which compromise DM in 
our universe today. Since ^ and ^ carry an asymmetry, it is natural (but not necessary) for 
their decay products to generate the baryon asymmetry, through 

— )• quarks , ^' — )■ x + quarks . (1-3) 

In this way, DM and baryons can arise from the late decay of a common parent particle, 
leaving a baryon asymmetry but symmetric DM. 

This mechanism of DM assimilation is general and can be implemented in many singlet 
DM models. In this work, we focus on the SUSY bino case for concreteness. We consider 
a minimal extension of the MSSM, where we only add a heavy vector-like chiral multiplet 
^/^^. The only new parameters with respect to the MSSM alone are the superpotential 
mass and soft terms for these new states. The new i?-odd degrees of freedom destroy and 
"assimilate" bino DM in the early universe through the process in Eq. (1.1). As long as 
^/ijrc have the right asymmetry and the correct quantum numbers, their late time decays 
can be responsible for generating the baryon asymmetry. We find that this mechanism works 
best when ^/^^ are at the TeV scale, implying that the LHC is able to probe this physics 
by looking for quasi-stable charged particles. 

^ Quasi-stability of and 4' requires that the triangle inequahty among the particle masses < m^+m* 
is satisfied (as well as its crossed versions). In the limiting case where 5' and ^I" are absolutely stable, the 
reaction in Eq. (1.1) is an example of semi-annihilation [41]. 
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{SU{3)c,SU{2)l,U{1)y) 


U{1)b 




(1,1, +1) 
(1,1,-1) 


+1 
-1 



Table 1: Additional superfield content compared to the MSSM for a simple model of DM 
assimilation. Here, we consider multiplets that only carry hypercharge, with generalizations 
shown in Table 3. 





Type 


Mass 


{SU{3)c,SU{2)l,U{1)y) 


U{1)b 


i?-parity 




Complex Scalar 


+ 2M 


(1,1,+1) 


+1 






Complex Scalar 


+ 2M 


(1,1,+1) 


+1 






Dirac Fermion 


M 


(1,1,+1) 


+1 


+ 



Table 2: Component fields and their properties for the heavy states. The masses are given 
in the limit 6 <^ M, with 6 defined in Eq. (2.8). For the relic abundance calculations, we will 
ignore the 6a/ mass contribution for simplicity. 

The remainder of this paper is structured as follows. In Sec. 2, we introduce a simple 
model of DM assimilation and mention possible generalizations. In Sec. 3, we study bino 
assimilation in the early universe, showing both semi-analytical and numerical results, and 
leaving calculational details to the appendices. We discuss the experimental implications for 
DM searches and collider measurements in Sec. 4, and conclude in Sec. 5. 

2 A Simple Model of Dark Matter Assimilation 

In this section, we present a simple model of DM assimilation involving new quasi-stable 
particles. These new particles carry a matter/antimatter asymmetry which we identify with 
SM baryon number, and their late time decay generates both DM and baryons. In Sec. 2.4, 
we describe a few generalizations of this model. 

2.1 Field Content and Lagrangian 

In our model, the MSSM is augmented by vector-like superfields ^/^'^ that carry hyper- 
charge and baryon number. The quantum numbers of ^/^^ are shown in Table 1, with 
alternatives discussed in Sec. 2.4. The leading superpotential terms consistent with the sym- 
metries are 

= Wmssm + M**'= . (2.1) 

We will work in the limit of a (nearly) pure bino, so the only MSSM multiplet relevant for 
our discussion is the hypercharge gauge multiplet, containing the hypercharge gauge boson 
= cos6w7^i + sinOwZfj, with coupling ay = 5'y/(47r), and the bino x with soft mass 
m^. We assume iZ-parity conservation throughout this paper, and tildes refer to i?-odd 
fields. This superpotential has a continuous U{1)b symmetry, under which ^ and have 
charges ±1, and we identify U{1)b with SM baryon number through decay operators shown 
in Sec. 2.3. 
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The physical degrees of freedom are summarized in Table 2. Including the relevant 
SUSY-breaking soft terms, the scalar masses are 



2 



Scalar ^ "(M^ + m^^^fj [ + V" " V'V'' + h.C. , (2.2) 



where ip (ip^) is the scalar component of ^ (^'^)- In the presence of the bM term, the scalar 
mass eigenstates are 

^1 = ^^, ^2 = ^7|-> (2-3) 

with masses 



mi,2 = + ± bM ■ (2.4) 

For simplicity, we will neglect the 6jvf terms for the relic abundance calculation in Sec. 3, but 
they are important when we discuss the scalar decays in Sec. 2.3. We will use the notation 
^ to refer to both ipi and ij;2 

$^{V^i,^2}, ^^^{i^li^l}, (2.5) 

since they have the same properties up to the small mass difference. 

It is convenient to combine the Weyl fermions ip/ijj^ in ^/^^ into a single Dirac field 

with mass term 

^fcrmion ^ "M^^ . (2.7) 

Ignoring the 6m term, the fermion/scalar mass splitting from the scalar soft terms is 



S = mscalar " mfermion = + m^^^^ - M ~ , (2.8) 

and we will assume that 6 > 0, such that the i?-odd scalars ^ are heavier than the R- 
even fermions ^. Since we expect to be comparable to msoftj the triangle inequality 
< mq, + rriij, is generically satisfied and the decay ^ — >• is forbidden. 

2.2 Matter /Antimatter Asymmetry 

With a mass M ~ few TeV, the abundance of \I' and ^ would be quickly depleted in the 
absence of some kind of primordial asymmetry. In order to have a large enough density of ^ 
and ^ to affect DM thermal freezeout, ^ and ^ will have a (shared) asymmetry dictated by 
SM baryon number U{1)b- In particular, we assume that the full SM asymmetry^ 



VB 



^ 8.6 X 10-11 (2.9) 



is stored in ^ and ^ in the early universe, and transferred to SM baryons via the late decay 
of ^ and ^. Here, (n-g) is the baryon (antibaryon) number density and s is the entropy 



^In contrast to the standard notation, we use a comoving definition of rjB- 
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Figure 1: Leading decay of the scalar ^' to a bino and SM quarks, through a squark. To 
avoid problems with the fermion ^ decays, we will always take the squark to be off-shell. 



density of relativistic degrees of freedom. We do not give an explicit model for ^'-ogenesis, 
but we assume that the asymmetry was established well before DM freezeout.'^ 

In Sec. 3, we will be interested in studying the comoving number densities Yi = rii/s. 
It is convenient to package together the scalar modes with the same baryon number 

= Yj +Y~ , Y-^ = Y-^ + Y-^ , (2.10) 

such that the total baryon asymmetry (before \I'/^ decay) is 

VB^Y^+Y^-Y^-Y^. (2.11) 

With this notation, there are no additional factors of two needed to account for the two 
complex scalars and one Dirac fermion. 

As the universe cools, antimatter will be efficiently destroyed, and during DM freezeout, 
r]B ^ Y^ +Y^. Since 5 in Eq. (2.8) is positive, Y^ < Y^ at all times (i.e. the i?-odd scalars 
are less abundant than the i?-even fermions). 

2.3 Late Time Decays 

The heavy ^ and ^ states will decay to the MSSM at late times. Since the scalars ^ are 
i?-odd, their decays necessarily terminate with the LSP, which in our scenario is the bino. 
These decays are responsible for bino DM production, yielding 

Ydm ^ 1$ , (2.12) 

similar to the super- WIMP scenario [42, 43]. Corrections to this relation arise because im- 
parity only ensures that ^ (^) decays will yield an odd (even) number of binos. Since bino 
DM requires Ydm ^ we must ensure that the i2-even states ^ do not produce very many 
bino pairs, which can happen by arranging the appropriate kinematics. We also assume that 
these decays do not violate U{1)b, such that at late times tjb is the true baryon asymmetry. 

With the quantum numbers of ^/'J''^, the leading dimension-5 decay operator involving 
MSSM fields is 

Wdecay = — T • (2.13) 

A,- 



decay 



^It is possible to enhance the effect of DM assimilation by assuming 'if carries a larger asymmetry, which 
cancels against an onti-asymmetry of SM baryons. 
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(a) (b) 

Figure 2: Leading decay channels for the fermion (a) three-body decay to SM quarks 
through a gluino loop; (b) five-body decay to bino pairs and SM quarks through tree-level 
(off-shell) squarks. As long as the squarks are kinematically inaccessible, the first process 
dominates the total width, as needed to forbid bino overproduction through ^ decays. 



Here, and are the electroweak singlet quark multiplets, and we are agnostic as to their 
flavor structure.^ Since the decay ^' — )■ "^x is kinematically forbidden, the dominant decay 
channel for the scalars ^ is 

■^^WWc^X, (2.14) 

through the Feynman diagram shown in Fig. 1. Note that the fejv/ mixing term in Eq. (2.2) 
is needed to ensure that both ipi and ^p2 have allowed decay modes. 
The leading decays of the fermion ^ are 

'^^iflFdP, (2.15) 
^i^uPWdPxX, (2.16) 

through the Feynman diagrams shown in Fig. 2. If the squarks were kinematically accessible, 
then decays like ^ ^ u'^u'^ dP ^ u'^ XX would dominate, leading to vast overproduction 
of binos. Thus, we must assume the MSSM squarks are not kinematically accessible (or at 
least phase space suppressed). That said, even if the squarks are off-shell, the process in 
Fig. 2b is still sizable, and we have to make sure that binos originating from this process 
are a subdominant component of the DM observed today. From the asymmetry we know 
that ~ ~ 10"^*^, and from observations we know the bino comoving density today is 
1~ ~ 10~^^, therefore the ratio between the partial widths for the two processes in Fig. 2 has 
to be less than 10~^. This ratio can be estimated by assuming all the final state particles 
are massless,^ and we assume that the squarks and gluinos have a mass comparable to M 
for simplicity. We find 

Jj^^ffi^ ~ 1^ ~2 X 10-^ (2.17) 
r — 72 (1-2 



*Since Eq. (2.13) involves an SU{3)c epsilon tensor, the two u"^ fields must have difi'erent fiavors. 
^Since the final state binos are not massless, Eq. (2.17) is somewhat of an underestimate. 
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where we have used the dimensionless volume of n-body phase space 



167r2 (n-l)(n-2) 



$(2) = 1. 

Svr ■ 



(2.18) 



Hence, as long as the squarks cannot be on-shell and the gluino is not too heavy, the loop- 



suppressed three-body decay ^ 



' dP can dominate over the dangerous five-body decay. 



In order to have the desired cosmology, ^ and ^ have to be long lived compared to 
DM freezeout in order to achieve DM assimilation, but short lived compared to Big Bang 
Nucleosynthesis (BBN) to ensure that the baryon asymmetry is correctly transferred to SM 
baryons. Again assuming that the squarks and gluinos have a mass comparable to M, we 
estimate the lifetime of these states as 



1 



ay 



2^A2_^ 245767r5 



7 X 10"^ sec 



1 TeV\' 



M 



-^decay \ 

108 TeV J 



a" 



M2 m2 



2M (167r2)2 A2^^^y 256^3 



2 X 10"' sec 



1 TeVy ( A. 

M 



decay 



/ V 108 TeV 



(2.19) 
(2.20) 



where the 1/I67r2 factors correspond to loop- or phase-space suppressions. For the right 
choice of Ajecay, these lifetimes are indeed long compared to DM freezeout (tf ~ 10"^ sec), 
but short compared to the start of Big Bang Nucleosynthesis (tsBN — 1 sec). Thus, we can 
treat ^ and ^ as stable for the DM freezeout calculation in Sec. 3, but still have standard 
BBN. It might be interesting to study modifications to BBN if A^ccay is taken to be larger. 

Finally, we note a bound on Eq. (2.13) coming from the fact that this operator can 
transfer the asymmetry from to SM baryons via high temperature scattering processes 
(as opposed to just decays). This would deplete the asymmetry in and make DM 

assimilation less effective.^ Requiring that the Hubble expansion rate is faster than the 
scattering rate T gives a bound on the reheat temperature 



Tr < 1.8 X 10' 



A2 

3 decay 

Mpi 



(2.21) 



Requiring > M to make sure that have a standard thermal history yields 



A, 



decay 



> 2.6 X 10^ TeV 



M 



1 TeV 



1/2 



(2.22) 



which is consistent with the story above. A similar constraint would also need to be satisfied 
in a full theory of ^'-ogenesis, as well as in the generalizations below. 

2.4 Generalizations 

There are a number of generalization of the model in Sec. 2.1 using the same superpotential 
in Eq. (2.1). As long as ^/^"^ carries baryon number (not necessarily ±1), has non-zero 
hypercharge, and has an allowed late-time decay to quarks, then we have a potentially viable 
model of (pure) bino assimilation. The -R-odd states must be heavier than the i?-even states 
to avoid overproduction of DM, and the i?-even state must have a suppressed decay to two 
binos, perhaps because of kinematics as in Sec. 2.3. A number of possibilities for ^/^^ are 
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{SU{3)„SU{2)l,U{1)y) 


U{1)b 


Hierarchy? 


Decay Operator 


(1,1, +1) 


+1 


^scalar ^ ^fcrmion 




(3,1, +4/3) 


+2/3 


^scalar ^ ^fcrmion 




(3,1, +2/3) 


+1/3 


^scalar ^ ^fcrmion 




(3,2,-1/6) 


-1/3 


^scalar ^ '^fcrmion 




(3,3,-1/3) 


-2/3 


'^scalar ^ '^fermion 




(1,4,-1/2) 


-1 


^scalar ^ ^fcrmion 


^qqq 



Table 3: A non-exhaustive list of alternative quantum numbers for the ^ multiplet in 
Sec. 2.1. The first line is the example studied in this paper. In each case, has the conjugate 
quantum numbers to ^ . The mass hierarchy is chosen to make sure that the i?-odd state is 
heavier than the i?-even state. The decay operators involve the SM quark multiplets, namely 
the electroweak doublet q and electroweak singlets u'^ /df^^ with unspecified flavor structure. 
The bilinear decay operators indicate small mixings with the MSSM states. 

shown in Table 3. In the cases where ^ /^'^ carry color or electroweak charges, there are 
additional diagrams compared to Sec. 3.1 involving the SU {2>)c and SU{2) l gauge multiplets. 

Beyond bino DM, there are other examples where singlet DM can be assimilated into 
a heavier particle in the early universe, and then produced from late time decays. Another 
SUSY example is a chiral superfield S which is a singlet under the SM gauge group, with 
the superpotential interaction with heavy fields 

Ws = Xs ^^^'^ + 7nsS^ . (2.23) 

The R-odd fermion singlino s DM would be assimilated into the -R-odd scalars ^. A key 
ingredient is an annihilation channel for the i?-odd states that conserves baryon number 

^ ^ (2.24) 

which can occur via i-channel gaugino and singlino exchange. The late-time decays which 
populate DM and baryons would occur through operators like Eq. (2.13). 

Finally, while the new states need to carry some kind of asymmetry for assimilation 
to be effective, it need not be the baryon asymmetry. Since the new states eventually need 
to decay, this asymmetry would persist until today, so the baryon asymmetry is the natural 
choice given phenomenological constraints.^ 

3 Bino Assimilation in the Early Universe 

In this section, we study in detail how bino assimilation happens in the early universe. After 
describing the relevant processes, we derive the Boltzmann equations for the species evolution, 
and present semi-analytical and numerical results. Various details appear in the appendices. 

®See however, footnote 3. 

^One amusing possibility would be if the new states could carry a DM asymmetry (as in asymmetric DM 
[40]). Assimilation could then be used to help eliminate the symmetric component. 
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Figure 3: Assimilation of the bino x, converting fermions ^ to scalars ^. The diagrams are 
(a) s-channel scalar exchange; and (b) t-channel fermion exchange. The field B denotes the 
hypercharge gauge boson. 

3.1 Relevant Processes 

We start the thermal history from early times, when the temperature of the universe is higher 
than the superpotential mass M but the baryon asymmetry is already stored in ^/^^. At 
this stage, the asymmetry in Eq. (2.11) is largely irrelevant, and particles and antiparticles 
democratically populate the universe. Once the temperature drops below the mass M, the 
antiparticles ^ and ^ quickly disappear as a consequence of annihilation, as shown in App. A. 
The antiparticles can be neglected at lower temperatures, leaving only \I' and \I' carrying the 
baryon asymmetry: 

rjB = Y^ + Y^. (3.1) 

Given the absence of antiparticles, there are only three classes of processes relevant for 
DM freezeout. The first relevant process is assimilation shown in Fig. 3, 

Assimilation: ^B, (3.2) 

where B is the hypercharge gauge boson.® Since the ^ and ^' states decay much later than 
DM freezeout and they satisfy the triangle inequality < + m^, this is an example 
of a semi-annihilation reaction [41]. In this process, the bino x g^ts destroyed, but its odd 
i?-parity is stored in the final state ^. In other words, DM has been assimilated. 
The second relevant process is destruction of i?-odd particles shown in Fig. 4, 

Destruction: x$ ^B. (3.3) 

This is also a semi-annihilation process, and the associated diagrams can be obtained by 
crossing symmetry from Fig. 3. The common features of these two process, shared with all 
semi-annihilation reactions, is that they are never phase space suppressed. Thus, they are 
very effective in the early universe at keeping the bino in thermal equilibrium, causing near 
total bino depletion. 

While Eq. (3.3) is effective at keeping the x iii thermal equilibrium, it is not the most 
effective way to eliminate i?-odd particles, since both assimilation and destruction are con- 
trolled by the same bino abundance. The third and last relevant process is i?-parity conversion 

^Though B is a linear combination of the mass eigenstates 7 and Z, we ignore the Z mass for simphcity. 
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Figure 4: Destruction of iZ-odd particles for: (a) s-channel fermion exchange; and (b) 
t-channel scalar exchange. 




Figure 5: Conversion of R-odd particles to R-even particles through t-channel bino exchange. 



shown in Fig. 5, 



Conversion : ^ ^ 



(3.4) 



As long as ^ are heavier than ^, this is most effective way to get rid of i?-odd particles 
in the early universe, despite the fact that this process is phase-space suppressed. In the 
semi-analytical study in Sec. 3.3, this conversion process will be solely responsible for the 
freezeout of the scalars. Finally, when the scalars and fermions decay at late times, the binos 
are repopulated with density Y^. 

3.2 Boltzmann Equations 

The Boltzmann equations describe the evolution of the comoving number densities of the 
bino X ™d the scalars ^. During DM freezeout, the comoving number density of fermions 
^ is given by baryon number conservation 



= r]B 



(3.5) 



so we need not have a separate Boltzmann equation for Y"^. 

We introduce a dimensionless time variable x and cross sections A,- 



^ A- 



s{T = M) 
H{T = M) 



rcl/i 



(3.6) 



- 10 - 



where T is temperature, H is the Hubble parameter, and s is the entropy density of relativistic 
degrees of freedom. The equihbrium comoving density for the bino is 



X 



eq 



9x 45 /m^ 





—X 




IM J 



(3.7) 



where g^^ = 2 and ii'2[a^] is the modified Bessel function. To express the equihbrium distribu- 
tion for the new degrees of freedom ^ and it is convenient to introduce the ratio between 
the equihbrium distributions which is not affected by the asymmetry 



(x) 



M 



using the parameter 5 from Eq. (2.8). The equihbrium distributions are 

1 



r^'^ix) 



^?(^) l + 4^(x) 



1 + r^'*(x) 

The Boltzmann equations, derived in detail in App. B, are 



dx 



ycqycq 
yeq -'ij, 



A 



dest 

x"^ 



Y~Y~ 



y^qycq 
^eq ^"I- 



dY^ 

dx 



y-eqy-^q 

^q— - ^X^* 



Y~ 



Ad 



est 



Y~Y~ 



y-eqyeq 
X I- 



Y^ 



Yl 




(3.8) 



(3.9) 



(3.10) 



(3.11) 



The effective cross section parameters can be expressed in terms of dimensionless cross sec- 
tions \i for the fundamental processes in our model 



A 



assim 
Adcst 



A. 



Aj$ ^ 



9 \ 



(3.12) 



Here, we are relying on the simplified scalar notation in Eqs. (2.5) and (2.10), and the 
resulting factors of two are tracked in detail in App. B. 



3.3 Semi- Analytical Solution 

In general, the two coupled Boltzmann equations in Eqs. (3.10) and (3.11) have no closed 
form solution, and we will solve them numerically in the next subsection. However, there is a 
particular region of parameter space where the number density evolution can be understood 
semi-analytically. This occurs if the assimilation and destruction processes are effective 
enough to keep x i^i thermal equilibrium, such that ^ and ^ freezeout when there are very 
few X around. In this case, we can focus on the ^ and ^ system alone, and the bino DM 
density today is determined directly from the freezeout density of the scalars ^. 
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If we can neglect the bino number density Y^, the assimilation and destruction processes 
are subdominant, and the number density of scalars ^ is approximately described by the 
conversion process alone 



Yi-K'^fiVB-Y^y . (3.13) 



dx X^ 

The above equation can be solved semi-analytically, in analogy with the Lee- Weinberg calcu- 
lation [6, 44]. At early times, the conversion reactions keep the ^ particles in thermal equilib- 
rium. Eventually, a freezeout point is reached where the reaction rate Tconv = f^~^(o"frci)conv 
cannot compete with the Hubble expansion rate K. We define the freezeout point x/- by 

if(x/) = crconv(x/) , (3.14) 

where c is a constant of order 1. This yields the relation cY^*^ Aconv = ^/j or equivalently 

CT]B Ac 



M 

Xf = Y log 



3M, 



(3.15) 



After freezeout, the factors in the Boltzmann equation are exponentially suppressed, and 
the evolution is straightforward. Ignoring the freezeout boundary condition, the asymptotic 
value for is 

^$(oo) = ^ . (3.16) 

Since the scalars ^ will eventually decay to the bino x, and since (almost) all the DM 
density comes from this decay, the mass density of DM today is 

Px,o = "^x'5o^$(oo) • (3.17) 
As usual we can express the DM density as a function of the critical energy density 



1.07xlO»QeV- 

^/^Mpi (cTWrcOconv M 



where the effective cross section {crvj-d) com is related to the parameter Aconv by Eqs. (3.6) and 
(3.12). This is the final semi-analytical result. The important difference between Eq. (3.18) 
and the familiar Lee- Weinberg result [6, 44] is the ratio m^/M. This arises because the 
particle freezing out is not the same as the one comprising DM today, analogous to the 
super- WIMP case [42, 43). 

To get a sense for the range of parameters that lead to the desired DM density, we 
compute the conversion cross section for our particular model in App. C. In the 6 <^ Ad limit 

27ra^ msoft 

(O-Urcl)conv ^ o 1^ ' {3.19} 

X 

where m^oft is related to 5 in Eq. (2.8). Combining Eqs. (3.18) and (3.19), we find 



where we have used = 106.75 and Mpi = 1.22 x 10^^ GeV. This is only an approximate 
result, valid in the 5 <^ AI limit and assuming that Y^ can be neglected. For reasonable 
values of and msoft, we obtain the desired relic density. The dependence on M comes in 
through Xf in Eq. (3.15). 
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Benchmark 


M 




"^soft 


A 


1.4 TeV 


180 GeV 


200 GeV 


B 


500 GeV 


390 GeV 


200 GeV 


C 


900 GeV 


180 GeV 


350 GeV 


D 


600 GeV 


390 GeV 


350 GeV 



Table 4: Four benchmark points that yield the observed DM rehc abundance JIdm^^ — 0.1. 
3.4 Numerical Results 

We now solve the Boltzmann system in Eqs. (3.10) and (3.11) numerically, using the thermally 
averaged cross section in App. C. We will first study four representative benchmark points 
given in Table 4. In all four cases, we have selected values for the scalar soft mass msoftj 
the superpotential mass M, and the bino mass m^^ that reproduce the observed relic density, 
f^DM^i^ — 0.1. For benchmarks A and C, the bino is much lighter than 'f/^, and for 
benchmarks B and D, the spectrum is more compressed. In Fig. 6, we show the evolution of 
the comoving number densities. 

Benchmarks A and B have a common value mgoft — 200 GeV, and their evolution is 
shown in Figs. 6a and 6b. For these two benchmarks, the bino distribution follows very closely 
its equilibrium value, which implies the semi-analytical analysis in Sec. 3.3 will be valid. The 
bino is kept in thermal equilibrium due to the very effective assimilation /destruction reactions 

^ $S , x$ ^ ^5 . (3.21) 

The above processes almost completely wash out the binos, which are "assimilated" into the 
scalars ^. Turning to the heavy states, if thermal equilibrium were maintained between ^ 
and ^, all the scalars \I' would eventually be suppressed, leaving no i2-odd particles around 
today, and thus no DM. However, the expansion of the universe prevents that from happening, 
and the scalar comoving density freezes out. From the semi-analytical result in Eq. (3.15), 
we know that the freezeout point scales as j;j ~ with a logarithmic dependence on 

the other variables, and this dependence of xj can be seen by comparing the departure from 
equilibrium between Fig. 6a and Fig. 6b. 

We next consider benchmarks C and D in Figs. 6c and 6d, which have a larger value 
"isoft — 350 GeV. We immediately notice a substantial difference from the previous case, 
such that the semi-analytical solution cannot be trusted. Increasing the soft mass yields a 
dramatic decrease in the \I' thermal equilibrium distribution from Eq. (3.9), because of the 
larger mass splitting appearing in the exponential.^ As a result, the overall interaction rates 
are smaller than in the previous case, and there is an early departure from thermal equilib- 
rium ("freezeout") while the x density is orders of magnitude above the ^ density. The 
subsequent evolution of the scalar particles ^ involves two competing effects. Immediately 
after freezeout, the assimilation process starts to regenerate since there is a large density 
of X that drives the assimilation reaction. Eventually, the x density decreases sufficiently that 
the destruction and conversion processes can reduce the number of ^ until the asymptotic 
value is reached. While the DM dynamics is richer and less straightforward for large values 

®The cross sections \i appearing in the Boltzmann system in Eqs. (3.10) and (3.11) are not appreciably 
affected by the higher soft mass value. 
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A: m^ = 180GeV, M = 1.4 TeV, d = 14 GeV 



B: »i^ = 390GeV, M = 500 GeV, = 39 GeV 




Figure 6: Numerical results for the Boltzmann system in the four benchmark scenarios from 
Table 4. The top row has mgoft — 200 GeV while the bottom row has rrisoft — 350 GeV. The 
left column has m^^ <^ M while the right column has ~ M. In all cases the parameters 
are chosen to reproduce the observed value of the DM density, 0dm ~ 0.1. We plot the 
comoving number densities for the bino x (blue) , Dirac fermion ^ (green) and scalars ^ (red) . 
We also plot the equilibrium distributions, using the same colors but dashed lines. The light 
red band shows the semi-analytical solution for the asymptotic scalars comoving density from 
Eq. (3.20), sweeping the range c = [0.1,2]. As discussed in the text, for benchmarks C and 
D, the semi-analytical result is not valid since the bino density is large at freezeout. At late 
times, the bino number density is given by Is,. 



of the msoft; the overall picture is left unchanged: the scalars ^ assimilate the bino x, and 
late-time scalar decays produce the DM observed today. 

Once the quantum numbers of ^/^'^ are fixed, the DM density depends on only three 
parameters: the superpotential mass M, the bino mass m^, and the soft mass mgoft- To 
further explore the allowed parameter space, we identify the region in the three-dimensional 
parameter space (M, m^, mgoft) where the bino DM density is within the WMAP 95% CL 
region [45], namely 0.0975 < JIdm^^ ^ 0.1223. In what follows, we plot slices of this three- 
dimensional parameter space obtained by fixing one parameter at a constant value. 

In Figs. 7a and 7b, we fix the soft mass to a constant value, and we identify the region 
in the (M, m^) plane which gives the observed DM density. In Fig. 7a we consider the case 
rrisoft = 200 GeV as in benchmarks A and B. We immediately observe that a heavier bino 
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Figure 7: DM relic density, considering slices through the (M, m^, mgoft) parameter space. 
The shaded blue regions identify the WMAP 95% CL region, 0.0975 < 0dm /i^ < 0.1223, and 
the contour lines correspond to constant DM density. The red dots indicate the benchmark 
points in Table 4. Top row: relic density in the (M, rriy^) plane for: (a) msoft = 200 GeV; (b) 
"^soft = 350 GeV. Bottom row: relic density in the (M, mgoft) plane for: (a) = 180 GeV; 
(b) = 390 GeV. 



requires a lighter superpotential mass. This makes sense, since we require a smaller comoving 
relic density y^(oo) for a heavier bino to get the same mass density, and the semi-analytical 
solution shows that Y^{oo) increases as M does. In Fig. 7b, we consider a higher value 
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of the soft mass, mgoft = 350 GeV, as in benchmarks C and D. As we discussed above, a 
higher soft mass leads to ear Her freezeout, and thus a higher DM density. We checked that 
for rrisoft > 400 GeV, the DM produced is always overabundant for the particular model in 
Sec. 2.1, but this 400 GeV limit can be evaded for some of the generalizations in Sec. 2.4. 

In Figs. 7c and 7d, we consider slices of constant bino mass, = 180 GeV and 
ruy^ = 390 GeV, respectively. We explore the allowed region in the (M, m-soft) plane. One 
interesting feature is that for certain values of and M, there are two options for mgoft that 
reproduce the observed DM density. For lower soft masses, freezeout occurs when <^Y^, 
as for the semi-analytical case. For a larger soft mass, freezeout occurs when Y^^Y^, and 
we have competition between assimilation and destruction/conversion. 



4 Implications for Experiments 

From the point of view of DM direct and indirect detection experiments, our model has the 
same phenomenology as the MSSM with (nearly) pure bino DM. Crucially, the relic density 
is different from the naive thermal expectation, since with assimilation, the bino abundance 
can be adjusted to the observed value by choosing the mass M appropriately. We saw in 
Sec. 3.4 that M had to be around 1 TeV for assimilation to be effective, and this mass range 
implies spectacular collider signals. Long-lived charged particles can be observed at the LHC, 
and if they are stopped, they exhibit distinctive multi-body decays. 

4.1 Dark Matter Searches 

For both direct and indirect detection experiments, the absence of any wino/higgsino content 
implies that the expected signals for pure bino DM are quite weak. Generic SUSY neutralino 
elastic cross sections for scattering off target nuclei was computed in several references (see 
e.g. Refs. [46, 47]), for both spin-independent and spin-dependent contributions. The pri- 
mary feature of the pure bino case is the absence of t-channel Higgs or Z exchange, which 
dramatically reduces the cross sections. 

The remaining contributions are from s-channel squark exchange. If we assume small 
A-terms (which is natural for the first- and second-generation quarks), then there is no 
mixing between the left- and right-handed squarks, and the spin-independent cross section 
from squark exchange is highly suppressed, as can be seen from the limiting expression 
in Ref. [48]. Thus, in a generic part of parameter space, the dominant spin- independent 
contribution comes from the small higgsino component, which allows for t-channel Higgs 
exchange, with a cross section [48, 49] 



^^v/proton,SI ^ 1-76 X lO'^ pb ( ^ ) ( ITl^ ) / \ , (4.1) 



Aan/3V /150 GeV^ fl TeV^ ^ 



\ 5 J \ niho J 



where is the Higgsino mass parameter. This cross section is just below current experimental 
limits (~ 10~^ pb) for small values of the bino/higgsino mixing angle. There is a spin- 
dependent contribution from s-channel squark exchange, which is well-approximated by [48] 



1 TeV^ ^ 



t^x/proton,SD ^ 7 X 10"" pb ( ) , (4.2) 

but this is well below current limits (~ 10^^ pb). 
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LHC (7 TeV), CalcHEP 3.1, CTEQ 6L 

0.1 
0.01 



0.001 

o-(pb) 
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10"^ 
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M (GeV) 

Figure 8: Leading order cross section for pair production of fermions ^ and scalars ^ as 
a function of M. These cross sections are given for the 7 TeV LHC, and compared to a 
LI fb"^ CMS bound on quasi-stable staus [52, 53]. Note that the efficiency and acceptance 
for the CMS search are not exactly the same as for pair-production, but the bound is 
indicative of the LHC sensitivity. The scalar cross section is smaller than the fermion cross 
section due to threshold suppression in electroweak pair production. 




The situation with indirect detection is quite pessimistic, since bino pair annihilation 
in our galaxy today is strongly suppressed. The only available channel is p-wave suppressed 
annihilation to leptons and quarks, computed in Refs. [46, 49, 50]. In the pure bino limit, the 
annihilation is usually dominated by leptons in the final state, because among the sfermions, 
the right-handed slepton is typically the lightest (assuming sfermion universality at a high 
scale) and has the highest hypercharge. The annihilation cross section is [15] 



^4..;. <,.3.3xlO-%b (— j (^5^) ,(4.3) 



which is well below current sensitivities (~ 1 pb). 

4.2 Quasi-Stable Charged Particles at Colliders 

Unlike the feeble signals predicted in DM searches, the heavy states provide spectacular 
collider phenomenology. From the relic abundance calculation, we know that M ~ O(TeV), 
so both ^ and ^ can be pair produced at the LHC. Since they are stable on collider time 
scales, they would be detectable as heavy stable charged particles [51]. 

Both CMS [52, 53] and ATLAS [54] already place bounds on heavy stable charged 
particles from anomalous ionization energy loss and time-of-flight measurements. For the 
model in Sec. 2.1, where only carry hypercharge, the most relevant study is from CMS 
[52, 53], which bounds stable staus using the 1.1 fb"^ data set. In Fig. 8, we show the cross 
sections for the fermions ^ and scalars ^ as a function of the mass M, calculated at leading 
order using CalcHEP 3.1 [55, 56] and CTEQ 6L parton distribution functions [57]. The CMS 
study assumed that the dominant production of staus came through SUSY cascade decays, 
whereas in our scenario, can only be pair produced through electroweak processes. This 
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affects somewhat the efficiency and acceptance for the search, but with that caveat in mind, 
M > 300 GeV is ahowed at present. 

The generahzations in Sec. 2.4 are potentially more interesting, especially for colored 
^/^^ which form i?-hadrons [58]. Because of the larger QCD pair production cross section, 
the LHC reach is extended to higher values of M. That said, the value of M typically has to 
be larger for colored ^/^^, since there are extra diagrams involving gluons/gluinos which 
contribute to the assimilation, destruction, and conversion processes. 

Most of the time, heavy stable charged particles simply exit the detector, but if they are 
produced with a small enough f3 factor, they can loose sufficient energy through ionization 
that they stop inside the detector [59, 60]. Searches for stopped exotics are currently being 
performed at the LHC [61, 62]. Moreover, a recent study in Ref. [63] suggests that ATLAS 
and CMS can measure the kinematics of the decay products. In the model in Sec. 2.1, 
^ has a three-body decay to jets while ^ has a four-body decay with three jets and one DM 
particle. It would be interesting to understand the extent to which these two decays can be 
distinguished. The generalizations in Sec. 2.4 would have different decay patterns, but the 
general expectation is two quasi-stable charged states that are nearly degenerate but have 
different decay rates, with the lighter state decaying to n SM particles, and the heavier state 
decaying to n SM particles plus one invisible DM particle. 

Finally, beyond the phenomenology of the heavy states, this scenario shares the same 
phenomenological features as the MSSM with the lightest neutralino as a pure bino. In 
general, robustly identifying the neutralino admixture is a challenging task [64], but the 
absence of a nearby chargino state would give some evidence that DM is a nearly pure gauge 
singlet (and thus must have non-standard cosmology). 

5 Conclusion 

Upcoming results from ground- and satellite-based DM experiments as well as results from 
the LHC have the potential to reveal important information about the dark sector. Assuming 
a positive detection of DM, one would ideally want to test whether or not DM is a thermal 
relic as expected from the WIMP paradigm. For that purpose, it is important to know 
whether the properties of DM measured today can be extrapolated to the early universe [65]. 

In this paper, we have introduced a new process called "assimilation" , where quasi-stable 
particles endowed with an asymmetry can affect the abundance of DM. These quasi-stable 
particles destroy DM particles but absorb DM number, and their subsequent annihilation 
and decay leads to the desired DM relic density. This mechanism allows pure singlet DM 
to have the correct thermal relic abundance, and opens up new possibilities to achieve pure 
bino thermal DM in the MSSM. 

While it may be surprising that long-lived particles can have such a dramatic affect 
on DM, it is worth remembering that the DM number density is rather small. The asym- 
metry rjB — 8.6 X 10^^^ is larger than the DM comoving number density today, as long as 
moM ^ 5 GeV. Thus, even small asymmetries can have a large effect on DM if appropriate 
interactions are present. This model offers an interesting twist on the growing literature 
trying to connect DM to the baryon asymmetry [66-84]. Here, the baryon asymmetry and 
symmetric DM are generated by the late-time decay of a heavy particle. The idea of tem- 
porarily storing the baryon asymmetry in a quasi-stable particle may open new directions 
for baryogenesis (see also Refs. [71, 80, 84]). 
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Despite the fact that DM is a singlet, this scenario does leaves an imprint on DM 
experiments. While pure bino DM has suppressed spin-independent scattering on nuclei, only 
a small higgsino component is needed to give a cross section within reach of direct detection 
experiments. The situation with indirect detection is less optimistic, since the same logic 
that leads to the expectation of bino overabundance also implies small annihilation rates in 
the DM halo. The most intriguing signature of this scenario comes at the LHC, which has 
the potential to observe TeV-scale quasi-stable charged particles. In this way, DM becomes 
visible through assimilation. 
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A Relic Density with an Asymmetry 

In this appendix, we show that in the presence of an asymmetry, antimatter particles are 
efficiently destroyed in the early universe. This justifies the fact that we neglect antimatter 
in the Boltzmann equations in Sec. 3. 

We study a simple toy example that generalizes the Lee- Weinberg scenario. We consider 
a massive particle ij: with mass and assume an asymmetry rj^ between particles and 
antiparticles 

^^^!!l^ = y^,_y^. (A.l) 

To study an example similar to our model, we consider the particles ip to be fermions charged 
under the hypercharge SM gauge group. The only possible channel for particle/antiparticle 
annihilation is 

ipip-^ BB , (A.2) 

where B is the hypercharge gauge boson. The Boltzmann equation describing the comoving 
number density evolution (see App. B for more details) is 



dY,i, X 



dx 



_ s{T = m^ 



where we have defined a dimensionless time variable x = m^/T. 

As usual, we start with equilibrium initial conditions at early times. Because particles 
and antiparticles have equal and opposite chemical potentials, the product Y^'^Y^'^ is not 
affected by the asymmetry, and it is equal to 

Y;''y;^ = Yl, , y^,o = ^^,^'K, [x] , (A.4) 

y*s 

where is the number of the internal degrees of freedom for the particle 'i/'- Using this 
condition, plus the asymmetry defined in Eq. (A.l), we find 
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1/^ = 8.56 X 10"", = 1.4 TeV 



t]^ = 8.56 X 10"", = 500 GeV 




Figure 9: Numerical results for the Lee- Weinberg scenario with an asymmetry, fixing rj^, = 
8.6 X 10~^^. The results for the number densities of ip (red) and ip (blue) are shown for two 
different particle mass: (a) = 1.4 TeV; (b) = 500 GeV. The equilibrium comoving 
densities are also shown in dashed lines. 



The Boltzmann equation in Eq. (A. 3) can be integrated numerically. The cross section 
for a ipip pair to annihilate to a pair of gauge bosons is [85] 

/ \ TT air , , r.\ 

We consider two different particle masses corresponding to the benchmarks A and B defined 
in Table 4, namely = 1.4 TeV and = 500 GeV, and motivated by the baryon 
asymmetry we choose t]^ ~ 8.6 x 10~^^. In Fig. 9, we plot the comoving number densities for 
both Ip and ip, as well as their equilibrium values. For small values of m^, or equivalently a 
larger annihilation cross section {o'v^ci)^^^bb^ antiparticle closely tracks the equilibrium 
distribution in Eq. (A. 4) and is negligible after x ~ 25. For larger value of m^, there is an 
early departure from equilibrium, but we see that even for m^, = 1.4 TeV the full picture 
survives: the antiparticles are destroyed efficiently after x ~ 25, and we can neglect them in 
our analysis. 



B Details of the Boltzmann System 

In this appendix, we derive the Boltzmann equations in Eqs. (3.10) and (3.11). We start 
from a quick review of the Boltzmann equations in an expanding universe.^'' We consider a 
generic massive particle V'a with mass iria, and we define a dimensionless time variable x, a 
comoving number density Y^, and dimensionless cross sections as follows 

_ma ^ _ "-a , _ s{T = mg) 

T s H{T = ma) 

where T is the temperature, s is the entropy density of the relativistic degrees of freedom, 
and H is the Hubble parameter. The quantity {(Jv^ci) ab^cd denotes the thermally averaged 



''A detailed derivation can be found in App. A of Rof. [41] 
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cross section, which is constant in the s-wave Umit. The Boltzmann equation for Ya is 



dx ^ ' ' 



A, 



ab^cd 



YaYb 



yeqyeq 



ab 



-^cq-^cq 
la If, 



YaYb - Y,Yd 



(B.2) 



where the sum over i runs over all possible reactions changing the number of V'a- The collision 
operator C" for a generic reaction ab — )• cd is given in two equivalent forms, either one of 
which is valid. It is typically more convenient to use the collision reaction with massless 
degrees of freedom in the final state such that the reaction is allowed at zero kinetic energy. 

From the arguments in App. A, we can safely ignore antiparticles in our scenario. That 
means that we have keep track of four species in our model: the bino Xi the complex scalars 
tpi and tjj2^ and the Dirac fermion ^ . Ignoring the hM term for simplicity, it is more convenient 
to work with ip and i^'^. These two complex scalars have the same mass, the same SM gauge 
quantum numbers, and the same interactions. Thus, if we assume the scalars have the same 
initial conditions, they will always have the same number density. For this reason we prefer 
to write everything in terms of the scalar comoving density Y^, defined as 



Y-f + Y^ 



Yj = Y^ 



(B.3) 



We now detail the possible reactions in our scenario. We start from the reactions 
changing the bino number density 



^B 



(B.4) 



The above reactions have the same cross sections when we replace tp"^ -H- ip, thus we have the 
following Boltzmann equation for the bino 



yfqycq 



dYx_sr^.^x_ o V^^B y^Y^ 5!^y- 

i L ip A L 

(B.5) 

Once we express everything in terms of = 2Y^, we recover the Boltzmann equation in 
Eq. (3.10). 

We now consider the scalar degrees of freedom y^. We first derive the Boltzmann 
equation for the scalar then we replace = 2Y^. The reactions changing the number of 

tjj are 



ipB ^ "IfX , i^X ^ ^ B , 'ipip ^ , 



(B.6) 



and we find the following Boltzmann equation for Y- 



dYj 
dx 



X 



x"^ 



yeqyeq 



yfq -'»/' 



Y-f-Y.Y^ 



X 



yfqycq 




(B.7) 
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With the replacement = 21^, we find the Boltzmann equation in Eq. (3.11). 

We finahy derive the Boltzmann equation for the Dirac fermion ^, which is not given 
in the main text. The reactions changing its number density are 



with corresponding Boltzmann equation 



-0X , 



dx 



A 



x'^ 



(B.9) 



As a cross check, we observe that the sum of the ^ and ^ Boltzmann equations gives 



d 



+ y<l. = const = r]B 



(B.IO) 



consistent with conservation of U{1)b- Hence, we only need to solve a system of two Boltz- 
mann equations, and we choose to solve for x and 



C Interactions and Thermal Cross Sections 

In this appendix, we give the interactions and thermal cross section for the model in Sec. 2.1. 
For simplicity, we consider the s-wave limit of the various interactions. 

The components of ^ and couple to the SM hypercharge gauge boson as 



where the covariant derivatives are 

D^ijj = d^il: + igyBf^tp 
Likewise, they couple to the bino via 



C D -\/2, 



V'^'^V''') x + x 



9Y 



Once we integrate out the auxiliary D fields, we also get a scalar potential 



1 



;9y 



\^\'' - |V'''P + MSSM fields 



(C.l) 

(C.2) 
(C.3) 

(C.4) 
(C.5) 



which is not relevant for our discussion. 

From the interactions above, it is straightforward to compute the thermally averaged 
cross sections (crui-ei) and use them to solve the Boltzmann equation system. We implement 
our model in FeynRules [86], using the superspace module described in Ref. [87]. We then 
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use the FeynRules interface to translate the model to CalcHEP format [55, 56], which we use 
to analytically compute the squared amplitudes \M.ab^cd\^ ■ In the (5 <C M limit, the squared 
matrix elements averaged over the initial and summed over the final polarizations are 



xi' - 



9'y 



Mmy{2M + my)2 



+ 



_5_ 

M 



2gY 7^. TTTTTT^ ' + O 



m^{M + m^){2M + m^Y 

(- 



_6_ 

M 



X 
/\/f3 



X 



(C.6) 



For a generic process ab — )• cd in the s-wave limit, we have 



11 



{an 



rel) ab—>-cd 



\M 



ab^cdl 



327r niami) 



+ 



mz. 



m 



2\2 



{ma + (rua + mt) 



1/2 



(C.7) 



The squared matrix element |A^a6->c(i| is function of the Mandelstam variables s and t, and 
for the s-wave amplitudes we have 



(ma + rrihY 



ma + mb 



rriamb 



(Ci 
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